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ABSTRACT 

In a previous paper Schaechter proposes using an extended Kalman filter to 
estimate adaptively the (slowly varying) frequencies and damping ratios of a 
large space structure. The present paper shovs that the time-varying gains for 
estimating the frequencies and damping ratios can be determined in closed-form 
so it is not necessary to integrate the matrix Riccati equations. After cer- 
tain approximations, the time-varying adaptive gain can be written as the. pro- 
duct of a constant matrix times a matrix derived from the components of the 
estimated state vector. This is an important savings of computer resources and 
allows the adaptive filter to be implemented with approximately the same effort 
as the non-adaptive filter. The success of this new approach for adaptive 
filtering has been demonstrated using synthetic data from a :wo mode system. 

I . INTRODUCTION 

Adaptive estimation and control techniques are being studied for their 
future application to the real-time control of large space structures, where 
uncertain or changing parameters may destabilize standard control system 
designs. In a recent paper Schaechter [l] proposes using an extended Kalman 
filter to estimate adaptively the (slowly varying) frequencies and damping 
ratios of a large space structure. For a system with N states and M (slowly var- 
ying) parameters the extended Kalman filter requires integration of an N+M by N-Ht 
nonlinear matrix Riccati equation to determine the covariance and gain for the 
filter. Schaechter introduces approximations which allow the integration of 
the nonlinear matrix Riccati equation to be replaced by integration of a smaller 
set of linear matrix equations. The N states of the system are estimated using 
constant gains determined off-line. The time-varying gains for estimating the 
(slowly varying) of M parameters are determined on-line by integrating an 
M by N set of linear matrix equations. 

The contribution of the work presented here is to show that the time- 
varying gains for estimating the (slowly varying) frequencies and damping 
ratios can be determined in closed-form so it is not necessary to integrate 
the M by N set of linear matrix equations. This is an important savings of 
computer resources and allows the adaptive filter to be implemented with 
approximately the same effort as the non-adaptive filter. In particular, 
after certain approximations the time-varying adaptive gain can be written 
as the product of a constant matrix times a matrix derived from the components 
of the estimated state vector. The constant matrix is determined off-line 
just as the constant gains for estimating the state are determined off-line. 
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The success of this new approach for adaptive filtering has been demon- 
strated on a computer simulation using synthetic data from a two mode system. 
Work in progress is applying the new approach to a much larger system using 
experimental data. The theoretical development and preliminary experimental 
results are presented in the paper. 


II. FORMULATION WITHOUT ADAPTIVE FILTERING 

The standard state variable formulation of the dynamic equations of motion 
are shown below where the dot indicates derivative, x is the state vector, u is 
the control vector, z is the measurement vector, and w and v are dynamic noise 
and measurement noise. [2] 


x ■ Fx + Gu + T w (1) 

z = Hx + v 

When the dynamic system is precisely known, a state estimator of the following 
form may be constructed where x indicates the estimate of the state x and K is 
the gain matrix. 

x * Fx + Gu +K(z - Hx) A (2) 

The differential equation for the estimation error x > x - x is obtained by 
subtracting Eq. (2) from Eq. (1). 

• 

x * (F - KH)x + Pw ~ Kv (3) 

The differential matrix equation for the covariance of the estimation error P 
follows where R and Q are from the covariance of tha measurement noise v an 
the dynamic noise v. 

P « E (%c T ) 

P - (F - KH)P + P(F - KH) T 

+ rqr T + krk t (4) 

The optimal gain matrix K is chosen to minimize the trace of the estimate 
error covariance to give the usual result 

K » PH 1 ^ 1 (5) 

Notice that for a precisely known dynamic system, the estimation gains may be 
precomputed, even in the event of a time varying system. The analysis used 
with the adaptive filter closely parallels the development without adaptive 
filtering. 


III. ADAPTIVE FORMULATION AND SOLUTION 

Adaptive control may be required when the model in Eq. (1) is unknown, 
uncertain, or dependent upon a changing system configuration. The modifica- 
tions that need to be made in Eq. (1) in order to include the effects of an 
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uncertain parameter are given below where the vector parameter a has a dynamics 
matrix C with dynamic noise w^, 

x ■ F(a)x + Gu + Tw 
a * Ca + w 

a 

z * Hx + v (6) 

As can be seen from Eq. (6), the system dynamics are now a function cf the vec- 
tor parameter a. In this formulation, the vector parameter a represents small 
changes from a nominal value co the average value of a is zero. These param- 
eters are assumed to be slowly varying so that they may be adjoined to the 
state vector. An adaptive state estimator may be written so both the state 
vector and the vector of parameters are updated using the measurements. 


x * F(a)x + Gu 4* K^(2 - Hx) 

a = Ca + K (z - Hx) 
a 


(7) 


Let the symmetric matrices P x and P a represent the covariance of th** *>rror in 
the estimates for x and a, respectively, and let the rectangular matrxx P ax 
represent the cross-covariance of the errors in the estimates of x and a. 

It is necessary to calculate these covariance matrices in order to determine 
the optimal gains l ^ and K a . The optimal gains are selected to minimize the 
trace of the covariance of the estimation error and have the following values. 


T -1 

K * P H X R 1 
x x 


f -I 

K * P H A R 
a 


( 8 ) 


Proceeding as before, and assuming the estimation error a * a - a is small, 
gives the vector differential equation for the error, 


x * (F - K H)x + 
x 


ts-1 


a + Tw - K^v 


(9) 


«-KHx + Ca + « - K v 

a a a 


The Ka^rix differencial equations for the covariance are: 
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where 

and 


dP x /dt - (F-KH) P x + P x (F-K x H) T 
+ Q + K RK T 

x XX 

+ < V > P a « +p .* I < p . i ) 1 

dP ix Mt - cp #x + p^CF-iyo 1 

+ P a (F a x) T - K a HP x 
T 

+ K RK 
a x 

dp /dt - CP + P C T + Q 
a a a a 

- \ W ax - P .x«.W T + *. “a 1 


F = 9F/8a 
a 

R, Q, and Q are covariances of v, Tw, and w 
a * 

(without delta function) . 


( 10 ) 


The remainder of this analysis will show approximations which can be used to 
reduce the computational effort needed to calculate the covariance matrices 
and the optimal gains when the covariance matrix P & is very sm*ll (of 
order t) and the covariance matrix P ax is also very small (of urder e) * The 
gain K a will be very small (of order e) because it is calculated from P ax , 

^he differential equation for the covariance matrix P x will involve somf small 
terms, but most of the terns are larger and constant. If the last two terms 
in the differential equation for P x are neglected (because they are small te:*ms 
of order e) , it is possible to calculate the steady-state constant value of the 
covariance P x . From the constant value of tha covariance P x the constant gain 
can be determined. As one might suspect, the constant gain has „he same 
value as it would have if there were no errors in estimating the parameters 

Because the covariance matrices P ax and P a are of order e, many of the tei 
In the differential equation for P a :?•* of order e squared. If the last 4 
terms in the differential equation for P a are neglected (because they are v 
small terms of order z squared), It is pi .sible to calculate the steady -state 
value of the covariance P fl (to order e). As one might suspect, the constant 
steady-state value obtained for P a is the same value which would have been 
obtained if K a were zero. 

All that remains is to calculate the time-varying covariance r d so .1 at the 
needed gain K a can be determined. Because the gain has been a\ to 

equal P X H T R"^, the last two terma £:• the differential equation foi D ay cancel 
out. For the remaining analysis it will be a turned there are N stale vari- 
ables so the first N/2 variables (designated he N/2 length * ^ctor x*) cor- 
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respond to mode position, an c the last N/2 variables (designated by the N/2 
length vect r x**) correspond to velocity of mode position. The differential 
equatioas c or the dynamics of the mode variables without any forcing or dis- 
turbing terms are presented below where A* corresponds to the damping terms 
(-25w) and A** corresponds to the frequency terras (- 0 )“) . Notice both A* and 
A** are diagonal N/2 by N/2 matrices. 


dx/dt - Fx 


dx*/dt = x** 


(ID 


d::**/dt - A*x*+A**x** 


Let there be N parameters in the vector a and arrange the order of the param- 
eters a so that the first N/2 parameters are the same as the elements of the 
dicgcnal matrix A* and the last N/2 parameters are the same as the elements 
of the diagonal matrix A**. Furthermore, assume the N-by-N symmetric co- 
variance matrix P a associated with these parameters is diagonal and composed 
of diagonal sub-matx* P a * and P a **. With these assumptions, the partial 
derivative can be w ** in a particularly simple way where x* and x** 

represent diagonal ut_ l JlC^S with the diagonal elements equal to the vectors 
x* and x** 


Fx 


F x 
a 


V V 


0 I 


x* ' 

A*A** 

«- J 


X** 



- P *** [x] 

a 


( 12 ) 
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One further assumption is that the dynamics matric C (for the parameters a) is 
diagonal and equal to the scalar Cq times the identity matrix I. With those 
assumptions, the differential equation for the cross covariance P ax can be 
written as fo^ows where x is a diagonal matrix made up of the elements of x. 


dP /dt = P (F-K H+C) T 
ax ax x 


+ (P a ***[x] ) T 


where 


and 


C = 


I 


0 

p *** = 

3 P * 

a 



( 13 ) 


The remainder of the analysis will deal with the cross-covariance matrix P^ 
which is the transpose of the covariance matrix P ax . The differential equa- 
tion for the cross-covariance P xa can be written as follows: 


where 


dP /dt = F* 
xa 


P 

xa 


+ P a***[ ; ] 


F* = F - K H + C 
x 


( 14 ) 


The linear matrix differential equation for P xa has particularly desirable char- 
acteristics. All the terms in the differential equation are known constants 
(because the gain K x and the covariance P a are known and constant) except for 
driving terms due to estimates of the state x. If the approximation is made 
that the derivative of the forcing terms x is equal to the dynamics matrix F 
times x, then, except for transient terms, the solution to the linear matrix 
differential equation for P xa can be written in closed form as a linear com- 
'na^ion of the torcing terms x. This is similar to the result in elementary 
linear differential equations where the general solution is composed of the 
sum of the homogeneous solution due to the unforced differential iquation and 
the particular solution due to the forcing function. 

Because the forcing function [x"| is a diagonal matrix, the first element Is 
the forci \g term for the first column of the solution for the matrix P a ^, the 
second element X£ is the forcing term for the second column of the matrix P ax , 
and so on. Let Pi be a vector which represents the i-th column of the matrix 
P^x* T^ e l insar matrix-vector differential equation for the i-th column can 
be written as follows where P a i»is a scalar which Is the i-th element of the 
diagonal matrix P a and x^ ’s a scalar which is the i-th element of x and P| is 
is the i-th col^nn of the matrix P a *** which is all zeroes except for entries 
equal to the diagonal elements of P a 
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(15) 


dP./dt = F* P, 
1 i 


+ P.* x. 
1 i 


The solution for the vector P^ is assumed to be composed of the sum of two 
vectors. The first vector is the constant vector times the scalar xj 
(corresponding to the estimate of the position of the mode) and the second 
vector is the constant vector times the scalar x^ (corresponding to the 
estimate of the velocity of the appropriate mode). 


P. = E. x. + G.x. (16) 

1 1 J 1 K 

where for i <N/2 then j = i and k = i + N/2 
for i >N/2 then j = i-N/2 and K * i 

The derivative of the vector P^ can be calculated directly if it assumed the 
derivative of the vector x is equal to F x with A* and A? being scalars which 
represent the j-th element of the representative diagonal matrices which make 
up F. 


dP^/dt = E^dx^/dt + G^dx^/dt 


= + G^A*Xj + Gj Aj **Xj [ 


07) 


Substituting the expression for the assumed form of the vector P^ and the 
expression for the derivative of the vector P^ into the differential equation, 
gives the following equations where 6^ is a discrete delta function which is 
unity if i equals j and zero otherwise. 


G i A j**j + <E i + G i A J** ) *k 


U8) 

* F * Vj + p * G i*k * S lj P j *- j + ! ik P k* *k 

Collecting all terms which multiply the scalar x. gives one. veccor equation 
and collecting all terms which multiply the scalar x k gives a second vector 
equation. The j are two vector equations and two unknown vectors and Gi. 


G iV ■ F * E i " Vi* 

E 4 + G lA) ** - F* Gj + 6 lk P k * 


(19) 
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The expression for E^ obtained fiom the second equation is substituted into the 
first equations to give a single equation with the unknown vector Gj 


G.A.* = F* (F*G , - G.A.** + 6,, P, *) + 5..P.* 
1 ] i i 3 ik k i] j 


( 20 ) 


Since nj* and Aj** are both scalars, it is possible to solve directly for the 
unknown vector G^ where I is the identity. 

G. = (IA.*+F*A.** - F*F*) _1 (6 . .P . *+6 .. F*P *) (21) 

l J J ij j ik k 

In the same way, the expression for obtained from the first equation is sub- 
stituted into the second equation to give a single equation with the unknown 
vector E-^. 


A . *E. = (F*-IA. **) (F*E. + . . P *) + A.*6..P * (22) 

J i J i ij j jikk 

It is again possible to solve directly for E^. 

E x - (I*.. + F* Aj .. - F*F*) -1 [ 6 ij (F *- IA J ** )P j* +S lk A j *V ] < 2 « 

Thus the two unknown vector quantities and G± can be determined from known 

quantities so the covariance vector F± and the approximation for the covariance 

matrix P can be determined, 
ax 
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IV. SIMULATION RESULTS WITH TWO MODES 


The new, simplified adaptive formulation was first tested with a single 
mode system. After encouraging results were obtained with one mode, a two-mode 
system was examined. The two-mode system used in the simulation studies is 
shown in Figure 1. 



FIG. I TWO-MODE SYSTEM 


The system consists of two masses, M, three springs, K, and three viscous 
dampers, B. For this study, M-l, K=l, and B»0.10. Control forces may be 
applied to both masses, random external forces disturb both masses, and noisy 
measurements of the position of both masses are available. The measurements 
are used for estimating the state vector, and for estimating the parameter 
vector. The differential equations representing this system are: 


Mx^ + 2Bx^ - BX 2 + 2Kx^ - “ f^ + 


Mx^ + 2 Bx 2 - Bx^ + 2 KX 2 - “ ^2 + W 2 


z i = x i + v i 
z 2 * x 2 + v 2 


(24) 


The natural frequencies and damping ratios of this system are: 

* 1 Cj * 0.05 

b>2 - 1.732 ? 2 - 0.0869 

where the low frequency mode is the common mode motion of the two masses. The 
spectral densities of both the process and measurement disturbances (Q and R) 
are 0.0163. Two hundred position measurements of both the masses were made 
during a sixty second computer simulation. This sixty second duration was 
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selected to assure about ten oscillations of the lowest frequency mode. The 
sample rate was selected to give about ten samples per cycle of the highest 
frequency mode. The value of the correlation time constant for the parameters 
that are to be estimated was 250 seconds (so C Q is 1/250). This value is much 
larger than the time constants of the system. The selection of a "large" 
value is important in order to allow the adaptive filter to average values 
over several cycles of the system. The fol’ wing table gives a summary of 
the test cases. In each case, both the standard, non-linear extended Kalman 
filter, and the simplified extended Kalman filter described in this paper 
were run in older to make comparisons. In all of the cases studied thus 
far, these two cases were indistinguishable, except for a small, initial 
transient. This transient effect is attributed to beginning the standard 
extended Kalman filter covariance integration with values slightly different 
from the steady state values. 


CASE 1 

u) ^ unknown 

known 

no noise 


unknown 

5 2 known 


CASE 2 

unknown 

unknown 

no noise 

— . — «| 

unknown 

£ 2 unknown 


CASE 3 

ojj unknown 

unknown 

noise 


a), ur.xnown 

4. 

t 2 unknown 



The results are shown i the following figures and are discussed below. 
In Case I the starting estimates for the natural frequencies were chosen 
to be 10% in error with estimated to be 0.9 (rather than 1.0, and u 2 
estimated to be 1.559 rather than 1.732). The damping parameters were 
exact, and no noise was present in the system. The results for the estimate 
of w l (Fig. 2) show that the modal frequency is very readily identified from 
the measurements, inspite of the 10% initial error in the estimate. As the 
system response diminishes, less information is available for updating the 
parameters. Consequently, with no new information coming into the system, 
the parameter estimate begins to return to its nominal value (0.9) with the 
selected time constant of 250 sec. The estimate of u 2 behaves similarly. 

In Case II, the objective was the same as in Case I with the additional 
problem of simultaneously estimating the damping parameters. The initial 
estimates of the damping parameters were zero. The results of the poor ini- 
tial guess of the damping parameter are evident in Fig. 3. The estimate 
of the modal frequency tends to be lightly damped, but in all other aspects, 
the estimate of w i appears to have the same features that were present in 
Case I. As has been found in past studies fa the estimate of the dampii $ 
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parameter itself is quite poor. This is due to the fact that the position 
measuxement contains very little damping information. 

Case III is identical to Case II with the addition of both process and 
measurement disturbances. Surprisingly, this case yielded the best results, 
as is evident in Figures 4 and 5. The effects of the noise are clearly 
visible in the figures. However, in contrast with the previous two cases, 
the process noise continues to excite the system after the transient effect 
of the initial conditions have subsided. The result is that the measurements 
continue to provide information on the parameters for the duration of the 
simulation. Since the higher frequency mode is more heavily damped, and is 
less perturbed by the external disturbance, the improvement in the natural 
frequency estimate of mode two is not as dramatic. 

CONCLUSIONS 

This paper has developed approximations which allow dramatic reductions 
in the on-line computational requirements of the extended Kalman filter. 
Numerical simulations of this technique have validated the approach for two 
simple spring-mass systems. It was found that the full non-linear extended 
Kalman filter and the closed-form adaptive filter developed in this paper gave 
virtually identical results. Work is currently in progress to apply this 
approach to a much larger system using experimental data, rather than simulated 
data. 
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FIGURE 4 Estimate of (CASE III) 
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FIGURE 5 Estimate of w 2 (CASE III) 
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